Interplay between electron-phonon couplings and disorder strength on the transport 

properties of organic semiconductors 
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The combined effect of bulk and interface electron-phonon couplings on the transport properties 
is investigated in a model for organic semiconductors gated with polarizable dielectrics. While the 
bulk electron-phonon interaction affects the behavior of mobility in the coherent regime below room 
temperature, the interface coupling is dominant for the activated high T contribution of localized 
polarons. In order to improve the description of the transport properties, the presence of disorder is 
needed in addition to electron-phonon couplings. The effects of a weak disorder largely enhance the 
activation energies of mobility and induce the small polaron formation at lower values of electron- 
phonon couplings in the experimentally relevant window < T < 300K. The results are 
discussed in connection with experimental data of rubrene organic field-effect transistors. 



In the last years, part of condensed matter research has 
focused on the realization of field-effect devices with cor- 
related oxides and organic materials. [1] This has opened 
up the possibility of injecting charge carriers without 
changing the chemical composition of the system, so 
that exploring phase diagrams of complex materials with 
small amount of disorder. On the other hand, the de- 
vices work with solid or liquid dielectric gates, therefore 
drawbacks arise due to different effects: for example, the 
lattice mismatch between bulk and gate, the interaction 
of bulk electrons with gate degrees of freedom, and the 
dielectric breakdown for strong electric fields. Finally, 
another important factor is the presence of traps at the 
interface between bulk and gate. 

Single crystal organic field-effect transistors (OFETs) 
have represented a step forward in the field of plastic 
electronics. They are even more important since the in- 
trinsic properties of organic materials can be explored. 
For example, single-crystal OFETs show a charge mo- 
bility at least one order of magnitude larger than that 
of thin films. Q Among them, the most promising and 
studied are those based on oligoacenes, like rubrene. [3[ 

Scaling laws of the mobility as a function of the di- 
electric constant of solid [3, [If and liquid @ gates have 
been discovered pointing out that the nearby dielectric 
has a strong influence. Actually, if the difference be- 
tween the dielectric constant of the bulk and the gate is 
small, at temperatures close or higher than 100K, the 
mobility \i of these systems exhibits a power-law behav- 
ior (p ~ T~ 7 , with 7 ~ 2). Q This conduction process 
cannot be simply ascribed to band transport, and 
it arises from the interaction of charge carriers with bulk 
low frequency intermolecular modes. On the other 

hand, if the dielectric constant mismatch is high, an in- 
sulating behavior is found with much smaller values of 
mobility at room temperature. [B| It has been proposed 
that the injected charge carriers undergo a polaronic lo- 
calization due to the interaction with the polarizable gate 
dielectrics. [!, HH However, accurate calculations have 
shown that the hopping barriers of localized polarons are 



too small with respect to the predicted values. [12J, |13| 
For example, the hopping barrier of rubrene transistors 
with a Ta-iO^ oxide gate is calculated to be 12.4 meV, 
a value much lower than 55 meV found in previous esti- 
mates. [f| Higher hopping barriers would require unreal- 
istically large electron-phonon (el-ph) couplings. 

Recently, the interplay between charge carrier coupling 
with intermolecular vibrations in the bulk and the long- 
range interaction induced at the interface with the di- 
electric gate has been investigated. [13] These interac- 
tions stabilize a polaronic state, the bond polaron, at 
much lower strengths than expected from the polar in- 
teraction alone. However, this state still requires values 
of the bulk el-ph coupling larger than those typically es- 
timated. Moreover, the transport properties with com- 
bined el-ph couplings have not yet been calculated and 
compared with experiments. 

In this work, we investigate the combined effects of 
bulk and interface el-ph couplings on the transport prop- 
erties in a regime of realistic values of bulk el-ph interac- 
tion at finite temperatures. We show that the bulk cou- 
pling affects the behavior of mobility below room tem- 
perature enhancing the coherent contribution, but it is 
ineffective on the incoherent high T contribution domi- 
nated by the interface coupling. Furthermore, we include 
the effect of a weak disorder due to bulk and interface 
traps. Actually, the interplay between long-range el-ph 
interactions and disorder effects is a largely unexplored 
issue. The disorder effects are able to enhance the hop- 
ping barriers of the activated mobility and to drive the 
small polaron formation at lower values of el-ph interac- 
tions. We point out that disorder is a key factor to get 
agreement with experimental data in rubrene OFETs. 

Model We study a one-dimensional Hamiltonian model 
with coupling to bulk and interface vibrational modes 14 1 



H — H e i + H { gl lk + H e i-Bulk + H^ter + H el-Inter- (1) 

In Eq. |lj, the electronic part H e i is 

H c ,i = -t ^2 (cjcj+i + h - c -) + X! e * ni ' ( 2 ) 
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where t is the bare electron hopping between the nearest 
neighbors on the chain with lattice parameter a, c- and 
Cj are the charge carrier creation and annihilation oper- 
ators, respectively, relative to the site Ri, ej is a local 
energy whose fluctuations in the range [— W, W] simu- 
late disorder effects in the bulk and at the interface with 
gate, rii — c\ci is the density operator. The transfer hop- 
ping t is estimated to be among 80meV and 120meV. |8| 
Due to the presence of shallow traps, Q disorder is not 
overwhelming and it is distributed according to the flat 
probability function Q ({e^}). 
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where xi and pi are the oscillator displacement and mo- 
mentum of the bulk mode at the site i, respectively, m 
the oscillator mass, k the elastic constant and uiBuik = 
\Jkjm the mode frequency. The intermolecular modes 
are characterized by small energies {huiBuik — 5 — 10 
meV, with h Planck constant) in comparison with the 
transfer hopping. [1, Q 

In Eq. {T]), H e i^Buik represents the term similar to 
the Su-Schrieffer-Heeger [lj| (SSH) interaction for the 
coupling to intermolecular modes 



Hei-Buik = a^2(x i+1 - Xi) (cjc i+1 



h.c. 



(4) 



with a coupling constant. In the adiabatic regime for 
bulk modes (HuiBuik *C t), the dimensionless quantity 
Asuifc = a 2 /Akt fully provides the strength of the electron 
coupling to intermolecular modes. The typical values of 
^Buik are in the intermediate coupling regime (in. this 
work, we take XbuIU — 0.1 suitable for rubrene). [9l [Io| 

In Eq. ([TJ, Hunter is the Hamiltonian of free interface 
phonons 



H Inter = ^/nter^4 a 9> 



(5) 



where ujxnter is the frequency of optical modes, a' q and 
a q are creation and annihilation operators, respectively, 
relative to phonons with momentum q. 

In Eq. ([1]), H e i- i n ter is the Hamiltonian describing the 
electron coupling to interface vibrational modes 

H e l-Intra = ^ M q me iqRi (<*« + a -q) > ( 6 ) 



where M q is the interaction el-ph term 



gflUJinter iqR^ Rq 



41 



R n + i?,- 
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with g dimensionless coupling constant, L number of lat- 
tice sites, and Rq cut-off length of the order of the lattice 



spacing a. In order to quantify this coupling, we use the 
dimensionless quantity Xinter = J2 q M q /2hu!i n t e rt- In 
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this work, we take Rq = 0.5a and hioi nter = 0.5£. 

In the following, we will use units such that a = 1, 
H = 1, e = 1, and Boltzmann constant ks = 1. We 
will analyze systems in the thermodynamic limit (up to 
L = 32) measuring energies in units of t ~ 100 meV. 

Method The temperature range relevant for transport 
properties is such that uiBuik < T, therefore the dynamics 
of intermolecular bulk modes can be accurately treated as 
classical. [lj| The electron dynamics is influenced by the 
bulk oscillator deformations whose probability function 
P ({xi}) is self-consistently calculated through a Monte- 
Carlo approach. [l(| For one particle in a lattice of L 
sites, the probability function is very close to that of 
free harmonic oscillators. [Tj} At a fixed configuration of 
bulk displacements {xi}, Eq.Q shows that the electron 
hopping depend on the specific pair. 

At a fixed configuration of displacements {x{\ and local 
energies {e m }> the resulting inhomogeneous model can be 
diagonalized using the following basis states 



Ux ({xi}, {e m }) \n qi , ., n qL >\i>, 
where the unitary transformation Uj {{xi}, {e m }) 
Ux({ Xl },{e m }) = 



(8) 



exp 



TTg f* tt x rt> { £ ™}) n 3 eiqRl ( a « - a -i) 



(9) 



provides an el-ph treatment equivalent to a generalized 
Lang-Firsov approach. 18] For long-range el-ph interac- 
tions, this approach is very accurate in the intermediate 
to strong coupling regime. [l|l [2(| Through the varia- 
tional parameters f q ({x;}, {e m }), interface phonons are 
displaced from their equilibrium position to a distance 
proportional to the el-ph interaction. The parameters 
/ g ({ji},{e m }) can be self-consistently calculated or, 
without losing accuracy, expressed as f q ({xi}, {e m }) — 
a({xi},{e m }) (M q /uj In ter)/(b{{xi},{e m })cos(q) + 1), 
with a({xi},{e m }) and b({xi},{e m }) configuration 
dependent variational quantities. [21] In Eq. 
\n qi , .,n qL > denotes the phononic basis in the occupa- 
tion number representation, and \i >= c\\0 >, with |0 > 
vacuum state, is the basis state for the electron. 

After the diagonalization at fixed displacements {xi} 
and energies {e m }, the solution is obtained by minimizing 
the free energy. A central quantity of this work is the 
mobility as a function of the temperature. The mobility 
\x ({xi}, {e m }) is determined starting from the real part of 
the conductivity Re[o~{uj)] ({xi}, {e m }), taking the limit 
of zero frequency and dividing for the particle density 
1/L of the system. In the linear regime, the real part of 
the conductivity is derived from the Kubo formula [221 ] 

Re[a(u>)} ({xi}, {e m }) = 
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2u L J Jte^UHmO)), (10) 

where f3 = 1/T, A(t) operator in Heisenberg representa- 
tion, j is the current operator 

3 = -iJ2 Sc ]+s c 3 h iJ d x i}) > ( n ) 

hj t s = t — a5(xj+s — Xj) denotes the generalized 

hopping and 8 = — 1,1 indicates the nearest neighbors. 

The procedure for the calculation of the expectation 
value of an observable O is the following: first, one makes 
the average (O ({xi},{e m })} over the eigenstates and 
eigenvectors of the corresponding inhomogeneous model, 
then, over the distributions P {{xi}) and Q ({ej}) making 
the integral 

i 3 

(12) 

by means of the Monte-Carlo procedure. [10| 

Results First, we analyze the behavior of the mobil- 
ity as a function of the temperature in the absence of 
bulk el-ph coupling and disorder. In Fig. [TJ we report 
the mobility for different values of Xinter analyzing the 
intermediate to strong coupling regime. With increas- 
ing Xinter, the mobility decreases in the range of values 
smaller than unity. The mobility is expressed in natural 
units, that is in terms of /in = ea 2 /h ~ 7cm 2 /(V ■ s), 
taking a — 1°A. 8] Therefore, the order of magnitude of 
the mobility is in agreement with data in rubrene OFETs 
for 1 < Xinter < 2. [5] Below room temperature, the co- 
herent metallic behavior is present only for Xinter < 1.3. 
With increasing temperature, for 1 < A/, iter < 2, there 
is a hopping behavior with the activation energy A: 
/i oc exp (— A/T). Even for the unrealistically large cou- 
pling Xinter =2.1, the activation energy A = OAlt ~ 41 
meV is still smaller than that estimated in some rubrene 
OFETs. Q This implies that further interactions have 
to be included in order to describe more accurately the 
mobility. 

The next step is to analyze the role of bulk el-ph cou- 
pling. However, when one combines the bulk and inter- 
face el-ph couplings, the main effect is to improve the 
description of mobility in the coherent low temperature 
regime, leaving the high T incoherent regime essentially 
unaltered for realistic bulk parameters. In the left panel 
of Fig. [5J we report the mobility as a function of the 
temperature for Xinter = 1.3 for bulk coupling Xsuik = 
and Xsuik =0.1 (appropriate to rubrene). The phys- 
ical situation corresponding to Xinter, = has already 
been investigated in previous works. [9|, [1(3, As shown 
in Fig. [3J /i shows a tiny coherent band-like behavior 
at low temperatures, but, with increasing T, it goes to- 
wards the activated behavior where the bulk coupling is 
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FIG. 1. Mobility (in units of /in = ea 2 /Fi) as a function of the 
temperature for different interface el-ph couplings Xinter at 
Asuifc = 0. The quantity A is the polaron activation energy. 

not effective (in the right panel of Fig. [5] the polaron ac- 
tivation energies A are very close to those with A Bulk = 
shown in Fig. [T]). [T3| Actually, the mobility interpolates 
between the behaviors with only one coupling. At low 
temperature, the diffusive contribution is ascribed to the 
modulation of the electron kinetic energy due to the bulk 
modes with SSH interactions. 
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FIG. 2. Mobility (in units of fig = ea 2 /h) as a function of the 
temperature for different bulk el-ph couplings at Xinter ~ 1.3 
(Left panel) and for different interface el-ph couplings Xinter 
at X Bu ik =0.1 (Right panel). 

In the right panel of Fig. [2l the mobility as a func- 
tion of the temperature is shown for different Xinter at 
Xsuik =0.1. A tiny metallic coherent behavior is present 
in the relevant temperature range for all the couplings. 
This coherent contribution is weakened with increasing 
Xinter , but it does not disappear. This element is in con- 
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trast with experimental data which show a more or less 
marked insulating behavior from 150 K to 300 K. There- 
fore, the theoretical prediction of mobility is not accurate 
even if bulk and interface el-ph couplings are active. 
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At W = l.Ot, the coherent metallic term is strongly weak- 
ened, the transition temperature to localized polaron is 
decreased, and the activated regime is more stable. 

In summary, the interplay between el-ph couplings and 
disorder effects influences the transport properties of or- 
ganic semiconductors gated with polarizable dielectrics. 
The el-ph interactions for realistic values of model param- 
eters are not sufficient to provide a satisfying description 
of mobility, so that we show that a weak disorder is able 
to fulfill the discrepancy between measured and calcu- 
lated transport data. Disorder effects are able to reduce 
the coherent metallic behavior at low temperatures and 
to provide high temperature activation energies in agree- 
ment with data in rubrene organic field-effect transistors. 
The important point is that not only the interaction of 
bulk electrons with gate degrees of freedom, but also the 
presence of shallow traps in the bulk and at the interface 
is relevant for the properties of field-effect devices. 



FIG. 3. Mobility fi (in units of fio = ea 2 /h) as a function 
of the temperature for different disorder strengths W and 
Xb-uM = 0.1 at Xinter = 1.3 (Left panel) and Xinter = 1.1 
(Right panel). The quantity A is the polaron activation en- 
ergy. 



In order to explain the experimental data, it is nec- 
essary to include also disorder effects. Indeed, there is 
evidence of traps in the bulk and at the interface with 
gates. @] In the left panel of Fig. El we show the mobil- 
ity as a function of the temperature with increasing the 
strength of disorder for A/„ ter = 1.3 and Xsulk = 0.1- 
There are two main results. The first one is related to 
the suppression of the coherent metallic behavior with in- 
creasing W. The second one is the strong enhancement 
of the activation energy A up to 0.67i ~ 67 meV even 
for the small amount of disorder W — 0.7i. Furthermore, 
the decrease of the magnitude of the mobility is not so 
marked. Therefore, weak disorder effects are able to pro- 
vide a very accurate description of the mobility resulting 
as key quantities for the interpretation of experimental 
data. 

Another effect of disorder is to drive the small polaron 
formation at lower el-ph couplings. This phenomenon has 
already been studied in a single impurity model, H^, 24 1 
but, now, it is addressed in a more realistic situation 
where the disorder is distributed on all the lattice sites. 
In the right panel of Fig. G2 we show that the mobility 
as function of the temperature changes from metallic to 
insulating small polaron behavior at the lower interface 
el-ph coupling Xinter — 1-1 and weak disorder W = 0.5t. 
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